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Abstract 



Periodic orbit quantization requires an analytic continuation of non-con- 
vergent semiclassical trace formulae. We propose two different methods for 
semiclassical quantization. The first method is based upon the harmonic in- 
version of semiclassical recurrence functions. A band-limited periodic orbit 
signal is obtained by analytical frequency windowing of the periodic orbit sum. 
The frequencies of the periodic orbit signal are the semiclassical eigenvalues, 
and are determined by either linear predictor, Pade approximant, or signal 
diagonalization. The second method is based upon the direct application of 
the Pade approximant to the periodic orbit sum. The Pade approximant al- 
lows the resummation of the, typically exponentially, divergent periodic orbit 
terms. Both techniques do not depend on the existence of a symbolic dy- 
namics, and can be applied to bound as well as to open systems. Numerical 
results are presented for two different systems with chaotic and regular clas- 
sical dynamics, viz. the three-disk scattering system and the circle billiard. 
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I. INTRODUCTION 



Semiclassical periodic orbit quantization is a nontrivial problem for the reason that 
Gutzwiller's trace formula for chaotic systems and the Berry- Tabor formula for 
integrable systems do not usually converge in those regions where the physical eigenenergies 
or resonances are located. Various techniques have been developed to circumvent the conver- 
gence problem of periodic orbit theory. Examples are the cycle expansion technique Q , the 
Riemann-Siegel type formula and pseudo-orbit expansions |>J, surface of section techniques 
H, and a quantization rule based on a semiclassical approximation to the spectral staircase 
[7|. These specific techniques have proven to be very efficient for individual systems with 
special properties, e.g., the cycle expansion for hyperbolic systems with an existing sym- 
bolic dynamics. The other methods mentioned have been used for the calculation of bound 
spectra of specific systems. 

Recently, an alternative method based upon filter-diagonalization (FD) has been intro- 
duced for the analytic continuation of the semiclassical trace formula |8|]9|]. The FD method 
requires knowledge of the periodic orbits up to a given maximum period (classical action), 
which depends on the mean density of states. The semiclassical eigenenergies or resonances 
are obtained by harmonic inversion of the periodic orbit recurrence signal. The FD method 
can be generally applied to both open and bound systems and has also proven to be a 



powerful tool, e.g., for the calculation of semiclassical transition matrix elements [ 10 1 and 



the quantization of systems with mixed regular-chaotic phase space [ITTJ. For a review on 
periodic orbit quantization by harmonic inversion see ||12||. 



In this paper we present two different techniques to make semiclassical periodic orbit sums 
convergent. The first method is an advanced version of harmonic inversion adapted to the 
special structure of periodic orbit signals given as sums of 5 functions The semiclassical 



signal, in action or time, corresponds to a "spectrum" or response in the frequency domain 
that is composed of a huge, in principle infinite, number of frequencies. To extract these 
frequencies and their corresponding amplitudes is a nontrivial task. In previous work p|,p|JT2 



the periodic orbit signal has been harmonically inverted by means of FD fl4Hl6| which is 
designed for the analysis of time signals given on an equidistant grid. The periodic orbit 
recurrence signal is represented as a sum over usually unevenly spaced 5 functions. A smooth 
signal, from which evenly spaced values can be read off, is obtained by a convolution of this 
sum with, e.g., a narrow Gaussian function. The disadvantages of this approach are twofold. 
Firstly, FD acts on this signal more or less like a "black box" and, as such, does not lend 
itself to a detailed understanding of semiclassical periodic orbit quantization. Secondly, the 
smoothed semiclassical signal usually consists of a huge number of data points. The handling 
of such large data sets, together with the smoothing, may lead to significant numerical errors 
in results for the semiclassical eigenenergies and resonances. Here, we propose alternative 
techniques for the harmonic inversion of the periodic orbit recurrence signal that avoid these 
problems. In a first step we create a shortened signal which is constructed from the original 
signal and designed to be correct only in a window, i.e., a short frequency range of the total 
band width. Because the original signal is given as a periodic orbit sum of 5 functions, this 
"filtering" can be performed analytically resulting in a band-limited periodic orbit signal 
with a relatively small number of equidistant grid points. In a second step the frequencies 
and amplitudes of the band-limited signal are determined from a set of nonlinear equations. 
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To solve the nonlinear system, we introduce three different processing methods, viz. linear 
predictor (LP), Pade approximant (PA), and signal diagonalization (SD). It is important to 
note that these processing methods would not have yielded numerically stable solutions if the 
signal had not first been band-limited by the windowing (filtering) procedure. Furthermore, 
this separation of the harmonic inversion procedure into various steps may elucidate a clearer 
picture of the periodic orbit quantization method itself, and even provides more accurate 
results than previous calculations P,[T2|] using FD. 

The second method is the direct application of the Pade approximant to slowly con- 
vergent and/or divergent periodic orbit sums |T7|] . In the former or the latter case, the PA 
either significantly increases the convergence rate, or analytically continues the exponentially 
divergent series, respectively. The PA is especially robust for resumming diverging series 



in many applications in mathematics and theoretical physics |18 |. An important example 
is the summation of the divergent Rayleigh-Schrodinger quantum mechanical perturbation 
series, e.g., for atoms in electric [W\ and magnetic |30| fields. In periodic orbit theory the 



PA has been applied to cycle-expanded Euler products and dynamical zeta functions |]2T] 



It should be noted that the Pade approximant is applied in both methods in a completely 
different context, namely, in the first method as a tool for signal processing ||l"3| , |22|| , and in 
the second for the direct summation of the periodic orbit terms in the semiclassical trace 
formulae. 

In Sec. m we present our first method to make semiclassical trace formulae convergent. 
After briefly reviewing the general idea of periodic orbit quantization by harmonic inversion 
in Sec. |II A| we construct, in Sec. [II B| , the band-limited periodic orbit signal which is ana- 
lyzed, in Sec. |11 C| > with the help of either LP, PA, or SD. In Sec. [TTl]we introduce our second 
method for semiclassical quantization, viz. the direct application of the Pade approximant 
to the periodic orbit sum. In Sec. [TV] we present and compare results for the three-disk 
repeller and the circle billiard as physical examples of relevance. A few concluding remarks 
are given in Sec. |V[ 



II. HARMONIC INVERSION OF PERIODIC ORBIT SIGNALS 

A. General remarks 

In order to understand what follows, a brief recapitulation of the basic ideas of periodic 
orbit quantization by harmonic inversion may be useful. For further details see |12| . 



Following Gutzwiller jE],|| one can write the semiclassical response function for chaotic 
systems in the form 

g sc (E)=g s c (E)+J2Aoe iSp ° , (1) 

po 

where gQ C (E) is a smooth function and S po and A po are the classical actions and weights 
(including phase information given by the Maslov index) of the periodic orbit (po) contri- 
butions. Equation ([!]) is also valid for integrable systems when the periodic orbit quantities 
are calculated not with Gutzwiller's trace formula, but with the Berry- Tabor formula |3| 
for periodic orbits on rational tori. The eigenenergies and resonances are the poles of the 
response function. Unfortunately, the semiclassical approximation ([!]) does not converge in 



3 



the region of the poles, and hence one is faced with the problem of the analytic continuation 
of g sc (E) to this region. 

As in previous work [P,PUI2"|| , we will make the (weak) assumption that the classical 
system has a scaling property i.e., the shape of periodic orbits is assumed not to depend on 
a scaling parameter, w, and the classical action scales as 

S P o = ws po . (2) 

In scaling systems, the fluctuating part of the semiclassical response function, 

g sc (w)=J2Aoe iwSp ° , (3) 

po 

can be Fourier transformed readily to yield the semiclassical trace of the propagator 

1 f+oo 

C sc (s) = — I g sc (w)e- lsw dw = ]T A po 5 (s - s po ) . (4) 



00 po 



The signal C sc (s) has 5 spikes at the positions of the classical periods (scaled actions) s = s po 
of periodic orbits and with peak heights (recurrence strengths) A po , i.e., C sc (s) is Gutzwiller's 
periodic orbit recurrence function. Consider now the quantum mechanical counterparts of 
g sc (w) and C sc (s) taken as the sums over the poles w k of the Green's function, 

s qm M = E — ^— - > (5) 

V w - w k + le 



1 r+oo 

C qm {s) = — / g^(w)e- lsw dw = -i£ d k e~^ s , (6) 

Z7T J-oo ^ 



oo 



with d k being the residues associated with the eigenvalues. In the case under study, i.e., 
density of states spectra, the dk are the multiplicities of eigenvalues and are equal to 1 for 
non-degenerate states. Semiclassical eigenenergies Wk and residues d k can now, in principle, 
be obtained by adjusting the semiclassical signal, Eq. (f|), to the functional form of the 
quantum signal, Eq. (H), with the {wk,dk} being free, in general complex, frequencies and 
amplitudes. This scheme is known as "harmonic inversion". The numerical procedure of 
harmonic inversion is a nontrivial task, especially if the number of frequencies in the signal 
is large (e.g., more than a thousand), or even infinite as is usually the case for periodic 
orbit quantization. Note that the conventional way to perform the spectral analysis, i.e., 
the Fourier transform of Eq. (f|) will bring us back to analyzing the non-convergent response 
function g sc (w) in Eq. ([!]). The periodic orbit signal @ can be harmonically inverted by 



application of FD fl4Hl6|l, which allows one to calculate a finite and relatively small set of 
frequencies and amplitudes in a given frequency window. The usual implementation of FD 
requires knowledge of the signal on an equidistant grid. The signal (§) is not a continuous 
function. However, a smooth signal can be obtained by a convolution of C sc (s) with, e.g., a 
Gaussian function, 

C(^) = ^EV^ V)>2 - (7) 
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As can easily be seen, the convolution results in a damping of the amplitudes, dk — > d k = 
c^exp (— wla 2 /2). The width o of the Gaussian function should be chosen sufficiently small 
to avoid an overly strong damping of amplitudes. To properly sample each Gaussian a 
dense grid with steps As « a/3 is required. Therefore, the signal (|7|) analyzed by FD 
usually consists of a large number of data points. The numerical treatment of this large 
data set may suffer from rounding errors and loss of accuracy. Additionally, the "black box" 
type procedure of harmonic inversion by FD, which intertwines windowing and processing, 
does not provide any opportunity to gain a deeper understanding of semiclassical periodic 
orbit quantization. It is therefore desirable to separate the harmonic inversion procedure 
into two sequential steps: Firstly, the filtering procedure that does not require smoothing 
and, secondly, a procedure for extracting the frequencies and amplitudes. In Sec. |11B| we 
will construct, by analytic filtering, a band-limited signal which consists of a relatively small 
number of frequencies. In Sec. [II C| we will present methods to extract the frequencies and 
amplitudes of such band-limited signals. 



B. Construction of band-limited signals by analytical filtering 

In general, a frequency filter can be applied to a given signal by application of the Fourier 
transform |22|-^J. The signal is transformed to the frequency domain, e.g., by application 



of the fast Fourier transform (FFT) method. The transformed signal is multiplied with 
a frequency filter function f(w) localized around a central frequency, Wq. The frequency 
filter f(w) can be rather general, typical examples are a rectangular window or a Gaussian 
function. The filtered signal is then shifted by — Wq and transformed back to the time domain 
by a second application of FFT. The filtered signal consists of a significantly reduced set 
of frequencies, and therefore a reduced set of grid points is sufficient for the analysis of the 
signal. This technique is known as "beam spacing" [^3|] or "decimation" of signals. 

The special form of the periodic orbit signal (BJ) as a sum of 5 functions allows for an 
even simpler procedure, viz. analytical filtering. In the following we will apply a rectangular 
filter, i.e., f(w) = 1 for frequencies w G [wq — Aw,wo + Aw], and f(w) = outside the 
window. The generalization to other types of frequency filters is straightforward. Starting 
from the semiclassical response function (spectrum) g sc (w) in Eq. (|3]), which is itself a Fourier 
transform of the "signal" (01), and using a rectangular window we obtain, after evaluating 
the "second" Fourier transform, the band-limited (bl) periodic orbit signal, 

1 rWo+Aw . 

CEOO = — / g sc (w)e- ls ^Mw 

Z7T Jwo-Aw 
1 rw +Aw . 

= — ]T A po / e lsW0 -< s - s *° )w dw 

Z7C p Q Jwq—Aw 

sin [(s - s po ) Aw] ■ Q 

The introduction of u>o into the arguments of the exponential functions in (Q) causes a shift 
of frequencies by — wo in the frequency domain. Note that C^s) is a smooth function and 
can be easily evaluated on an arbitrary grid of points s n < s max provided the periodic orbit 
data are known for the set of orbits with classical action s po < s max . 
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Applying now the same filter as used for the semiclassical periodic orbit signal to the 
quantum one, we obtain the band-limited quantum signal 



1 rw +Aw . 

Cft m (s) = — / g (im {w)e- ls{w - Wo) dw 

27T J wo— Aw 



K 



-i Y, 4e" i( " o)s , \w k - w \ < Aw . (9) 



k=l 



In contrast to the signal C qm (s) in Eq. the band-limited quantum signal consists of a 
finite number of frequencies w^, k = 1, . . . , K, where in practical applications K can be of 
the order of ~ (50-200) for an appropriately chosen frequency window, Aw. The problem of 
adjusting the band-limited semiclassical signal in Eq. to its quantum mechanical analogue 
in Eq. (|^) can now be written as a set of 2K nonlinear equations 

K 

Cft(m-) = Cn = -i d k e~'< nr , n = 0, 1, . . . , 2K - 1 , (10) 

fe=i 

for the 2K unknown variables, viz. the shifted frequencies, w' k = Wk — w , and amplitudes, d k . 
The band-limited signal now becomes "short" as it can be evaluated on an equidistant grid, 
s = nr, with relatively large step width r = n/ Aw. It is important to note that the number 
of signal points c n in Eq. flT0p is usually much smaller than a reasonable discretization of the 
signal C^ c (s) in Eq. (0), which is the starting point for harmonic inversion by FD. Therefore, 
the discrete signal points c n = C^f(nT) are called the "band-limited" periodic orbit signal. 
Methods to solve the nonlinear system, Eq. (|10|) , are discussed in Sec. |I Q below. 

It should also be noted that the analytical filtering in Eq. @ is not restricted to periodic 
orbit signals, but can be applied, in general, to any signal given as a sum of 5 functions. An 



example is the high resolution analysis of quantum spectra |L2^g5| , g6| , where the density of 
states is g(E) = ^ n 8(E — E n ). 



C. Harmonic inversion of band-limited signals 

In this section we wish to solve the nonlinear set of equations 

K 

c n = Y,d k z n k , n = 0,l,...,2tf-l, (11) 

k=l 

where z k = exp (— iw' k r) and d k are, generally complex, variational parameters. For nota- 
tional simplicity we have absorbed the factor of — i on the right-hand side of Eq. ( [TO] ) into the 
rffe's with the understanding that this should be corrected for at the end of the calculation. 
We assume that the number of frequencies in the signal is relatively small (K ~ 50 to 200). 
Although the system of nonlinear equations is, in general, still ill-conditioned, frequency 
filtering reduces the number of signal points, and hence the number of equations. Several 
numerical techniques, that otherwise would be numerically unstable, can now be applied 
successfully. In the following we employ three different methods, viz. linear predictor (LP), 
Pade approximant (PA), and signal diagonalization (SD). 
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1. Linear Predictor 



The problem of solving Eq. (O) has already been addressed in the 18th century by 
Baron de Prony ||27|| , who converted the nonlinear set of equations (|TT|) to a linear algebra 
problem. Today this method is known as linear predictor (LP). Our method strictly applies 
the procedure of LP except with one essential difference; the original signal C sc (s) is replaced 
with its band-limited counterpart c„, = (^(nr). 

Equation ([11]) can be written in matrix form for the signal points c n+ i to c n+ K, 





( 


V c n+K ) 





v n+l 



y n+l 
Z K 



\ 



n+K 

Z K 



(dl \ 

\d K ) 



(12) 



From the matrix representation (|T2"D it follows that 



\z'{ 



n+K 



K 



y n+K 



V 1 


( Cn+1 \ 


J 


\ C n+K J 



K 



+k 



(13) 



k=l 



which means that every signal point c n can be "predicted" by a linear combination of the 
K subsequent points with a fixed set of coefficients a^, k = 1,...,K. The first step of 
the LP method is to calculate these coefficients. Writing Eq. ( |T3"D in matrix form with 
n = 0, . . . , K — 1, we obtain the coefficients at as the solution of the linear set of equations, 



( ci • • • c K \ 

\CK ■■■ C2K-1 J 



(ax \ 
\a K J 



( C ° ^ 



The second step consists in determining the parameters Zk in Eq. (|TT|). Using Eqs. 
(pTT|) we obtain 



K K K 

k=l 1=1 k=l 



n+k 



(14) 



and 



(15) 



and thus 



K 

E 

k=l 



K 



1=1 



z k z k 



4 = 



(16) 



Equation (1TB) is satisfied for arbitrary sets of amplitudes dk when Zk is a zero of the poly- 
nomial 



K 



. 



(17) 



The parameters z k = exp (— iw' k T) and thus the frequencies 
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log(*fc) 



(18) 



are therefore obtained by searching for the zeros of the polynomial in Eq. (|17D . Note that this 
is the only nonlinear step of the algorithm, and numerical routines for finding the roots of 
polynomials are well established. In the third and final step, the amplitudes dk are obtained 
from the linear set of equations 



E « 

k=i 



n 



0. 



K- 1 . 



(19) 



To summarize, the LP method reduces the nonlinear set of equations ((TTJ) for the varia- 
tional parameters {zk, dk} to two well-known problems, i.e., the solution of two linear sets 
of equations (|14]) and ( |T9| ) and the root search of a polynomial, Eq. (|T7|), which is a nonlinear 
but familiar problem. The matrices in Eqs. ([14]) and (|19"D are a Toeplitz and Vandermonde 
matrix, respectively, and special algorithms are known for the fast solution of such linear 
systems ||28|| . However, when the matrices are ill-conditioned, conventional LU decomposi- 
tion of the matrices is numerically more stable, and, furthermore, an iterative improvement 
of the solution can significantly reduce errors arising from numerical rounding. The roots of 
polynomials can be found, in principle, by application of Laguerre's method f28fl . However, 



it turns out that an alternative method, i.e., the diagonalization of the Hessenberg matrix 

/ 



A 



a K-l 


a-K-2 


. "1 


ap 


\ 


a K 


a K 


a K 


a K 


1 





■ 










1 


■ 








V 



o 



o 



(20) 



o J 



for which the characteristic polynomial P{z) = det[A — zl] — is given by Eq. flT7|) (with 
a = —1), is a numerically more robust technique for finding the roots of high degree 



(K > 60) polynomials [EH] 



2. Pade Approximant 

As an alternative method for solving the nonlinear system (|ll]) we now propose to apply 
the method of Pade approximants (PA) to our band-limited signal c n . Let us assume for 
the moment that the signal points c n are known up to infinity, n = 0, 1, . . . oo. Interpreting 
the c n 's as the coefficients of a Maclaurin series in the variable z~ l , we can then define the 
function g(z) = J2^=o c n z ~ n - With Eq. (]TT| ) and the sum rule for geometric series we obtain 

oo K oo K j p / \ 

g(z) = E = E 4 E ( = E - ttt\ • ^ 

n=0 k=l n=0 k=l Z Zk ^K{ z ) 

The right-hand side of Eq. ( ET]) is a rational function with polynomials of degree K in the 
numerator and denominator. Evidently, the parameters z& = exp (— iw' k r) are the poles of 
g(z), i.e., the zeros of the polynomial Qk{ z )- The parameters dk are calculated via the 
residues of the last two terms of (0). We obtain 



S 



ZkQ K [Zk) 

with the prime indicating the derivative d/dz. Of course, the assumption that the coeffi- 
cients c n are known up to infinity is not fulfilled and, therefore, the sum on the left-hand 
side of Eq. (|2"T|) cannot be evaluated in practice. However, the convergence of the sum 
can be accelerated by application of PA. Indeed, with PA, knowledge of 2K signal points 
Co, • • • , C2K-1 is sufficient for the calculation of the coefficients of the two polynomials 

K K 

P K {z) = ]T b k z k and Q K (z) = £ a k z k - 1 . (23) 

k=l k=l 

The coefficients a k , k = 1, . . . , K are obtained as solutions of the linear set of equations 

K 

C n = a k c n+k , Tl = 0, . . . , K - 1 , 
k=l 



which is identical to Eqs. (13) and (14) for LP. Once the a's are known, the coefficients b k 



are given by the explicit formula 

K-k 

bk= a k+m c m , k = 1, . . . , K . (24) 

It should be noted that the different derivations of LP and PA yield the same polynomial 
whose zeros are the z k parameters, i.e., the z k calculated with both methods exactly agree. 
However, LP and PA do differ in the way the amplitudes, d k , are calculated. It is also 
important to note that PA is applied here as a method for signal processing, i.e., in a 



different context to that in Sec. Ill, where the Pade approximant is used for the direct 



summation of the periodic orbit terms in semiclassical trace formulae. 



3. Signal Diagonalization 



In Refs. |ff4]JT6[1 it has been shown how the problem of solving the nonlinear set of equa- 
tions (|TTD can be recast in the form of the generalized eigenvalue problem, 

UB fe = z k SB k . (25) 

The elements of the K x K operator matrix U and overlap matrix S depend trivially upon 



the c n 's [16 



U, 



ij 



c i+j+l j 



Sij 



K- 1 . 



(26) 



Note that the operator matrix U is the same as in the linear system (|l"4D, i.e., the matrix 
form of Eq. ( |TB| ) of LP. The matrices U and S in Eq. ( |23| ) are complex symmetric (i.e., 
non-Hermitian), and the eigenvectors B k are orthogonal with respect to the overlap matrix 
S, 



f B k \S\B k i 



N k 5, 



kk' 



(27) 



9 



where the brackets define a complex symmetric inner product (a\b) = (b\a), i.e., no complex 
conjugation of either a or b. The overlap matrix S is not usually positive definite and therefore 
the iVys are, in general complex, normalization parameters. An eigenvector Bk cannot be 
normalized for Nk = 0. The amplitudes dk in Eq. (|Tl~D are obtained from the eigenvectors 
Bk via 



dk = — 



K-l 



n=0 



(28) 



The parameters z k in Eq. (JTT|) are given as the eigenvalues of the generalized eigenvalue 
problem (p5f) , and are simply related to the frequencies w' k in Eq. ( |10|) via z fc = exp(— iw' k T). 

The three methods introduced above (LP, PA and SD) look technically quite different. With 
LP the coefficients of the characteristic polynomial ([17]) and the amplitudes dk are obtained 
by solving two linear sets of equations (|14|) and (|1|). Note that the complete set of zeros Zk 
of Eq. (|i~7l ) is required to solve for the dk in Eq. (]l~9|). The PA method is even simpler, as only 
one linear system, Eq. (0), has to be solved to determine the coefficients of the rational 
function P K (z)/Q K (z). Finding the zeros of Eq. (|l7j) provides knowledge about selected 
parameters Zk, and allows one to calculate the corresponding amplitudes dk via Eq. (p2|). 
The SD method requires the most numerical effort, because the solution of the generalized 
eigenvalue problem (^5|) for both the eigenvalues Zk and eigenvectors Bk is needed. 

It is important to note that the three techniques, in spite of their different derivations, are 
mathematically equivalent and provide the same results for the parameters {zk,dk}, when 
the following two conditions are fulfilled: the nonlinear set of equations ( pTT]) has a unique 
solution, when, firstly, the matrices U and S in Eq. (^) have a non-vanishing determinant 
(det U 7^ 0, det S ^ 0), and, secondly, the parameters Zk are non-degenerate (zu ^ zy for 
k ytz k'). These conditions guarantee the existence of a unique solution of the linear equations 
(|I4|) and ([19]), the non-singularity of the generalized eigenvalue problem (p5|), and the non- 
vanishing of both the derivatives Q' K (zk) in Eq. (|22| ) and the normalization constants Nk 
in Eqs. ( |27| ) and (|28|). Equation ( |TT]) usually has no solution in the case of degenerate z k 
parameters, however, degeneracies can be handled with a generalization of the ansatz (|Tl~]) 
and modified equations for the calculation of the parameters. Here, we will not further 
discuss this special case. 

While the parameters z^ in Eq. ( |TT] ) are usually unique, the calculation of the frequencies 



w k via Eq. (|18|) is not unique, because of the multivalued property of the complex logarithm. 
To obtain the "correct" frequencies it is necessary to appropriately adjust the range Aw 
of the frequency filter and the step width r of the band-limited signal (]T0|). From our 
numerical experience we recommend the following procedure. The most convenient approach 
is to choose first the center wq of the frequency window and the number K of frequencies 
within that window. Note that K determines the dimension of the linear systems, and 
hence the degree of the polynomials which have to be handled numerically, and is therefore 
directly related to the computational effort required. Frequency windows are selected to be 
sufficiently narrow to yield values for the rank between K pa 50 and K pa 200. The step 
width for the band-limited signal should be chosen as 

(29) 



3 max 

i 



2K 
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with s max being the total length of the periodic orbit signal. The relation Zk = exp (— iw' k r) 
projects the frequency window w' G [—Aw, +Aw] onto the unit circle in the complex plane 
when the range of the frequency window is chosen as 

Aw = - = . (30) 

T "5 max 

When calculating the complex logarithm with arglogz G [ — 7r, +tt], Eq. flTBp provides the 
"correct" shifted frequencies w' k and thus the frequencies w k = w + w' k . 

To achieve convergence, the length s max of the periodic orbit signal must be sufficiently 
long to ensure that the number of semiclassical eigenvalues within the frequency window 
is less than K. As a consequence the harmonic inversion procedure usually provides not 
only the true semiclassical eigenvalues but also some spurious resonances. The spurious 
resonances are identified by low or near zero values of the corresponding amplitudes dk 
and can also be detected by analyzing the shifted band-limited signal, i.e., signal points 
ci, . . . , C2K instead of Co, . . . , o^k-i- The true frequencies usually agree to very high precision, 
while spurious frequencies show by orders of magnitude larger deviations. 

The harmonic inversion method introduced above will be applied in Sec. [IV A| to the 



periodic orbit quantization of the three-disk scattering system, and the semiclassical reso- 
nances will be compared with results obtained by the cycle-expansion technique ^,29,30 



and the direct application of the Pade approximant to periodic orbit sums discussed in the 
next Section |ITT[ 

III. SEMICLASSICAL QUANTIZATION BY PADE APPROXIMANTS TO 

PERIODIC ORBIT SUMS 

The method presented in the previous Section is based on signal processing (harmonic 
inversion) of the periodic orbit signal. In this Section we introduce our second method 
for semiclassical quantization, based on the direct application of the Pade approximant to 
periodic orbit sums. 

The PA to a complex function f(z) is defined as a ratio of two polynomials and can be 
computed from the coefficients a n of the Maclaurin expansion of f(z), i.e., a power series 



/(*) = E 



n=0 

with finite or even zero radius of convergence in z ||20|| . However, Eq. (|1|) does not have 
the functional form of a Maclaurin power series expansion of g(E) in the energy E. Even 
disregarding this limitation, a direct computation of the PA to the sum in Eq. ([I]) would be 
numerically unstable due to the typically large number of periodic orbit terms. Nevertheless, 
considering E as a parameter, g(E) can be rearranged and written as a formal power series 
in an auxiliary variable z, 

g(z;E)=Y: E A po {Ey s ^ E)/n z n ee £a n (£)^ , (31) 

n fj, po =n n 

where the maximum value of n required for convergence of the PA is relatively small com- 
pared with the number of periodic orbit terms. Of course, the arrangement (|31~D of the 
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periodic orbit sum as a power series is not unique, and expansions similar to Eq. (|3T|) can 
be used for other ordering parameters n of the orbits, e.g., the cycle length in systems with 
a symbolic code. However, if no symbolic dynamics exists, the sorting of orbits by their 
Maslov index is natural both physically and as a way to introduce an integer summation 
index and will be justified in Sec. [IV] by the successful numerical application of the method 
to a regular system without symbolic dynamics. Note that the PA to the periodic orbit sum 
cannot be applied without any ordering parameter, i.e., when no symbolic dynamics exists 
and all Maslov indices are zero, which is the case, e.g., for the Riemann zeta function as a 
mathematical model for periodic orbit quantization The true value g(E) of the periodic 
orbit sum is obtained by setting z — 1 in Eq. (|3l|), i.e., g{E) = g(l;E). In such a case, 
we have a point PA which is given as a ratio of two polynomials in z whose coefficients are 
non-polynomial functions of E all at a fixed value of z. The usual implementation of the PA 
as a ratio of two polynomials in z whose coefficients are computed, e.g., via the Longman 
algorithm [[H]] would be advantageous if g(z; E) were required for many values of z. In the 
present case, at each given energy E only one fixed value z = 1 is needed and the PA is 
most efficiently computed by means of the recursive Wynn e-algorithm ||32|| . 

To briefly describe the ^-algorithm, we introduce a sequence of partial sums {A n } which 
converges to (or diverges from) its limit A as n —>■ oo. In the case of divergence, A is called 
the 'anti-limit' of {A n } , as n — > oo. Further, let F be a transformation which maps {A n } 
into another sequence {B n }. The mapping F will represent an accelerator, i.e., the sequence 
{B n } will converge to the same limit A faster than {A n } if the condition 

Br, — B 

~T 1 

A n — A 

is fulfilled as n — > oo. In addition, the same F can be applied to wildly divergent sequences 
{^4 n }. Only nonlinear mappings can simultaneously accomplish both goals to accelerate 
slowly convergent and induce convergence into divergent sequences. The transformation F 
will be nonlinear if its coefficients depend on A n , e.g., the so-called e-algorithm of Shanks 
|33| , whose mapping F is the operator e k which converts the sequence {A n } into {B n } via 

e k (A n ) = B Kn = [n + k/k] (32) 

with n > and n > k. This is the well-known Aitken A 2 -iteration process (i.e., the simplest 
PA, [1/1]) extended to higher orders k. The general term in the /cth-order transform B^ n 
of A n can be computed efficiently from the stable and recursive e-algorithm of Wynn p2 
viz., 



e s (A m ) = 47 } = [m + s/s] , (33) 



where 



effi = e& 1) + (m+1) m J rn,s>0 (34) 
el - el 

with etj = 0, 4 m) = An, £ 2s+i = l/ e s(AA m ), and where A is the forward difference 
operator: Axj When {A„} is the sequence of partial sums of a power series, 
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the two-dimensional array £2™ yields the upper half of the well known Pade table [m/s]. 
However, the e-algorithm need not necessarily be limited to power series. 

The procedure to apply the above PA to semiclassical quantization by summation of 
periodic orbit terms works as follows. For a given system we calculate the periodic orbits up 
to a chosen maximum ordering parameter n < n max , where n can be but is not necessarily 
related to the Maslov indices of orbits. Note that this set of orbits usually differs from the set 
of orbits with classical action S po < S max , which is required for periodic orbit quantization 
by harmonic inversion in Sec. ||. From the periodic orbit amplitudes A po (including the 
phases exp(— i|/i po ) given by the Maslov indices) and actions S po , we compute the partial 
sums (e.g., with n < /i max the Maslov index): 

A n = £ A po {Ey s ^' h . (35) 

The sequence {A n } of partial sums is used as input to the e-algorithm, Eq. fl3"lp, to obtain 
a converged value g(E) of the periodic orbit sum ([J). The semiclassical eigenenergies or 
resonances are given as the poles of g(E) and are obtained by searching numerically for 
the zeros of the reciprocal function l/g(E). Such a search requires the evaluation of, e.g., 
Ap (E) and S po (E), at complex values of E. This can be done in a straightforward manner 
for the scaling systems considered in this paper. The number of the poles of g(E) is not 
constrained by the size of the sequence {A n } of partial sums since our PA is a ratio of two 
non-polynomial functions of E. 

It is important to note that the Pade approximant is applied here in a completely different 
context than in the first method where it is used as a tool for signal processing. In Sec. [I! 
we have shown that the three techniques LP, PA and SD used for the harmonic inversion of 
band-limited signals are mathematically equivalent. However, the two methods for periodic 
orbit quantization introduced in Sections and [TXT| are not necessarily equivalent, and it may 
well be that one or the other method is more appropriate for the semiclassical quantization 
of a given physical system. Numerical results for two physical systems with completely 
different dynamical properties will be presented in the next Sec. [TV|. 



IV. RESULTS AND DISCUSSION 



In this section we want to demonstrate the efficiency and accuracy of the methods in- 
troduced above by way of two examples: the three-disk repeller as an open physical system 
with classically chaotic dynamics and the circle billiard as a bound regular system. Both 
systems have previously been investigated by means of FD [ptp|, |i2| , p4| -|36|] . The three-disk 
repeller has also served as a prototype for the development and application of cycle expan- 
sion techniques P J29|J30[| , and for this system we will compare the convergence properties of 



three different methods, viz. harmonic inversion, Pade approximant and cycle expansion. 



A. The three-disk repeller 



As the first example we consider a billiard system consisting of three identical hard 
disks with radius R, displaced from each other by the same distance d. This simple, albeit 
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nontrivial, scattering system has served as a model for the development of the cycle expansion 
method |4, 21, 29, 30] and periodic orbit quantization by harmonic inversion |]8|j9|JT2]] . The 
three-disk scattering system is invariant under the symmetry operations of the group C^ v , 
i.e., three reflections at symmetry lines and two rotations by 27r/3 and 4ir/3. Resonances 
belong to one of the three irreducible subspaces A\, A 2) and E [[IT]. As in most previous 
work we concentrate on the resonances of the subspace A\ for the three-disk repeller with 
R = 1 and d = 6. In billiards, which are scaling systems, the shape of periodic orbits 
does not depend on the energy E, and the classical action is given by the length L of the 
orbit (S po = hkL po ), where k — |k| = y/2ME/h is the absolute value of the wave vector 
to be quantized. We have calculated all periodic orbits with Maslov index /j, po < 30, which 
corresponds to the set of orbits with cycle length n < 15. 



1. Harmonic inversion of the periodic orbit signal 

We first calculate the semiclassical resonances of the three-disk repeller by the method 
introduced in Sec. |T|, viz. harmonic inversion of the periodic orbit signal [see Eq. @] 

C SC (L) = J2Ao5(L-L po ) . 

po 

[Setting h — 1, we use s = L in what follows.] In Fig. |l]a we present the periodic orbit signal 
C SC (L) for the three-disk repeller in the region < L < L max = 35. The signal is given as a 
periodic orbit sum of delta functions S(L — L po ) weighted with the periodic orbit amplitudes 
A po - The groups with oscillating sign belong to periodic orbits with adjacent cycle lengths. 
Signals of this type have been analyzed (after convolution with a narrow Gaussian function, 
see Eq. ([7])) by FD in Refs. PHT2"]- We now illustrate the harmonic inversion of band-limited 
periodic orbit signals by LP, PA and SD. 

In a first step, a band-limited periodic orbit signal is constructed as described in Sec. 
I [T|. As an example we choose K = 100 as the rank of the nonlinear set of equations (|TU), 



and k = 200 as the center of the frequency window. The width of the frequency window is 
given by Ak = 27rK/L max = 2007r/35 ~ 18.0. The step width of the band-limited signal is 
r = AL = L max /2K = 0.175. The signal points c n = Cfi(L = nAL), with n = 0,...,2K 
are calculated with the help of Eq. (R) and presented in Fig. |T|b. The solid and dashed lines 
are the real and imaginary parts of Cy(nAL), respectively. The modulations with spacings 
7r/A/c result from the superposition of the sinc-like functions in Eq. (|8]). 

The band-limited periodic orbit signal C^(nAL) can now be analyzed, in a second step, 
with one of the harmonic inversion techniques introduced in Sec. [II Q , viz. LP, PA or SD. 
The resonances obtained by LP are presented as plus symbols in Fig. |Tjc. The dotted lines 
at Re k = 182 and Re k = 218 show the borders of the frequency window. The two symbols 
very close to the border on the left indicate spurious resonances. 

A long range spectrum can be obtained by choosing several values Wo for the center of the 
frequency window in such a way that the windows slightly overlap. By analyzing a periodic 
orbit signal similar to that in Fig. |l]b but with an increased signal length, L max = 55 we have 
calculated the semiclassical resonances of the three-disk repeller in the range < Re k < 250. 
It turns out that they are even more accurate than those obtained previously |P,|T2"|| using 



FD. For more details see Ref. |13 ]. Part of the resonances in the range 25 < Re k < 65 are 
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marked as squares in Fig. 0. The comparison of resonances in Fig. ^| obtained by various 
semiclassical quantization methods will be discussed below. 



2. Pade approximant to the periodic orbit sum 



We now apply our second method introduced in Sec. |TJ, the PA to the periodic orbit 
sum to the same system as discussed above, viz. the three-disk repeller with R = 1 and 
d — 6. We have calculated the partial sums {A n } of periodic orbit terms (see Eq. |35]) 
using all periodic orbits with cycle length n < 15, which corresponds to the set of orbits 
with Maslov index /i po < 30. The sequence of the partial sums {A n } of periodic orbit 
terms in Eq. (|35| ) converges for wave numbers k above the borderline Im k = —0.121557 
Q which separates the domain of absolute convergence of the periodic orbit sum from the 
domain where analytic continuation is necessary, but strongly diverges deep in the complex 
plane, where the resonance poles are located. This is illustrated in Fig. |3| for two different 
wave numbers k. The dashed line and the plus symbols in Fig. |3]a show the convergence 
of the sequence {A n } at k = 150 — O.li. What is plotted is the absolute value of the 
difference between two consecutive terms, e = \A n — A n ^\\. As can be seen this sequence 
is slowly convergent, and about two to three significant digits are obtained at n = 15. 
The convergence can be accelerated using the PA, as is seen by the solid line and squares in 
Fig. ||a showing the error values e = \A^ A — A^\\ for the sequence of the Pade approximants 
{A^ A } to the periodic orbit sum. The periodic orbit sum has converged to six significant 
digits already for n = 10. The situation is even much more dramatic in the deep complex 
plane (Im k < —0.122), e.g., at k — 150 — 0.5L Here, the sequence of the partial sums {A n } 
of periodic orbit terms exhibits exponential divergence, as can be seen from the dashed line 
and plus symbols in Fig. |3|b. Nevertheless, this sequence converges when subjected to the 
PA implemented through the ^-algorithm (solid line and squares in Fig. [3]b). 

The resonances of the three-disk scattering systems have been obtained by a numerical 
two-dimensional root search in the complex fc-plane for the zeros of the function l/g(k), 
where g(k) is the PA to the periodic orbit sum. The subset of the semiclassical resonances 
in the range 25 < Re k < 65 are marked by the plus symbols in Fig. |2|. They agree 
well with the squares obtained by harmonic inversion of the periodic orbit signal. For the 
resonances shown in Fig. |2] we will now discuss the applicability and limitations of the 
different semiclassical methods to some more extent. 



3. Harmonic inversion and Pade approximant vs. cycle expansion 

The three-disk scattering system discussed above has purely hyperbolic classical dynam- 
ics and has been used extensively as the prototype model for the cycle expansion techniques 
P J21|j29|j30|| . As has been shown by Voros [p^j , Gutzwiller's trace formula for unstable peri- 



odic orbits can be recast in the form of an infinite and non-convergent Euler product over 
all periodic orbits. When the periodic orbits obey a symbolic dynamics the semiclassical 
eigenenergies or resonances can be obtained as the zeros of the cycle expanded Gutzwiller- 
Voros zeta function. Unfortunately, the convergence of the cycle expansion is restricted, due 



to poles of the Gutzwiller- Voros zeta function [^TJ . The domain of analyticity of semiclassical 
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zeta functions can be extended |39|,|40[| resulting in the "quasiclassical zeta function" |f40|,|30 



which is an entire function for the three-disk repeller. This approach allows one to calculate 
semiclassical resonances in critical regions where the Gutzwiller-Voros zeta function does 
not converge, at the cost, however, of many extra spurious resonances and with the rate of 
convergence being slowed down tremendously |30 . 



It is very interesting and illustrative to compare the convergence properties of the various 
methods which have been applied for periodic orbit quantization of the three-disk repeller, 
viz. the two methods introduced in Sections [Q| and |T| of this paper and the cycle expansion 
technique. We will demonstrate that the harmonic inversion method and the PA to periodic 
orbit sums provide semiclassical resonances in energy regions where the cycle expansion of 
the Gutzwiller-Voros zeta function does not converge. With the limited numerical accuracy 
of harmonic inversion by FD applied in Ref. ||, the semiclassical resonances of the three-disk 
repeller in the region Im k < —0.6 were somewhat unreliable. The analysis of band-limited 
periodic orbit signals introduced in the present paper now allows us to calculate semiclassical 
resonances of much improved accuracy even deep in the complex plane. 

Fig. ^] presents the results of the three different methods in the region 25 < Re k < 65. 
The resonances calculated by harmonic inversion and the PA to periodic orbit sum are 
marked by the squares and plus symbols as already mentioned above. The crosses label the 
resonances obtained by the cycle expansion of the Gutzwiller-Voros zeta function |5C]. The 
dotted line in Fig. § indicates the borderline, Im k = —0.121557 Q, which separates the 
domain of absolute convergence of Gutzwiller's trace formula from the region where analytic 
continuation is necessary. For the two resonance bands slightly below this border the results 
of all three semiclassical quantization methods are in perfect agreement. The dashed line in 
Fig. |2| marks the borderline of absolute convergence of the Gutzwiller-Voros zeta function, 
at Im k = —0.699110 The Gutzwiller-Voros zeta function provides several spurious 



resonances which accumulate at Im k ~ —0.9, i.e., slightly below the borderline of absolute 
convergence (see the crosses in Fig. The resonances in the region Im k < —0.9, especially 
those belonging to the fourth band, are not described by the Gutzwiller-Voros zeta function 
but are obtained by both the harmonic inversion method and the Pade approximant to the 
periodic orbit sum (see the squares and plus symbols in Fig. 0, respectively). 



B. The circle billiard 

As the second example we now demonstrate the applicability of our semiclassical quan- 
tization methods to the circle billiard. This system has previously been investigated as a 
model for periodic orbit quantization of integrable systems |fH|]l2|j, however, to the best of 



our knowledge, has not yet been treated by the cycle expansion technique [Q or pseudo-orbit 
expansion p. 

The exact quantum mechanical eigenvalues E = h 2 k 2 /2M of the circle billiard are given 
by zeros of Bessel functions J\ m \(kR) = 0, where m = 0, ±1, ±2, , ... is the angular momen- 
tum quantum number and R is the radius of the circle. The semiclassical eigenvalues can 



be obtained by an Einstein-Brillouin-Keller (EBK) torus quantization |^3[ resulting in the 
quantization condition 

/ I TTl I f 3 \ 

kR\J 1 — [mj kR) 2 — |m| arccos = it ( n + - J (36) 
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where n = 0, 1,2, . . . is the radial quantum number. In the following we choose R = 1. 
States with angular momentum quantum number m^O are twofold degenerate. 

The periodic orbits of the circle billiard have the form of regular polygons. They can 
be labeled by two integer numbers m r and with the restriction m r > Im^ which can 
be shown to be identical with the number of sides of the corresponding polygon and its 
number of turns around the center of the circle, respectively ||41|| . Some examples are given 
in Fig. [|. The angle between two bounces is a rational multiple of 2tt, i.e., the periods L po 
are obtained from the condition 

L po = 2m r sin7, (37) 

with 7 = TTm^/rrir. Periodic orbits with m r ^ 2m^ can be traversed in two directions and 
thus have multiplicity 2. For the amplitudes A po of the circle billiard, the Berry- Tabor 
formula || for integrable systems yields 

•v=v!5 e " li(3m ' +4) ' (38) 

with /i po = 3m r being the Maslov index. 

The periodic orbit quantities can be used to set up the semiclassical recurrence signal 
(d) for the circle billiard which can then be analyzed by harmonic inversion to extract the 
eigenenergies. Detailed comparisons between results obtained by harmonic inversion and the 



EBK torus quantization fl3"6| ) are presented in Refs. [|i2| , |3"3| -|3"6|| and show excellent agreement. 
The harmonic inversion method can even be generalized, firstly, to the harmonic inversion of 
cross-correlated periodic orbit sums f35||H|| which allows to significantly reduce the number 



of orbits required for semiclassical quantization, and, secondly, to include higher order h 
corrections in the periodic orbit sum. We do not report these results here but refer the 
reader to the literature for details. In what follows we want to apply our second method 
introduced in Sec. |TTT| to the circle billiard. 

It follows from Eq. ( p8[) that the Pade approximant to the periodic orbit sum should be 
calculated with an ordering parameter n = m r — 1, i.e., 

A n {k) = J2 Aoe ikLp ° , n = 1, 2, . . . , (39) 

m r <n 



with the lengths L po and amplitudes A po given by Eqs. (|37|) and fl38|), respectively. We 
included all periodic orbits (m^m^) with m r < 100 in the calculation of the function g(k) 
which is the Pade approximant to the sequence A n (k). The real and imaginary parts of 
l/g(k) are presented as solid and dashed lines respectively in Fig. |5|. The zeros of the 
function l/g(k) agree perfectly to at least seven significant digits with the exact positions 
of the semiclassical eigenvalues obtained from Eq. (^) and marked by the squares in Fig. |5|. 



V. CONCLUSION 

We have introduced two different methods to make semiclassical periodic orbit sums con- 
vergent. The first method is based on the harmonic inversion of band-limited periodic orbit 
signals. The characteristic feature of this method is the strict separation of the two steps, 
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viz. firstly, the analytical filtering of the periodic orbit signal and, secondly, the numerical 
harmonic inversion of the band-limited signal by application of either linear predictor (LP), 
Pade approximant (PA), or signal diagonalization (SD). The separation of these two steps 
and the handling of small amounts of data compared to other "black box" type signal pro- 
cessing techniques provides a transparent approach to harmonic inversion as a semiclassical 
quantization method thus opening the possibility to an easier and deeper understanding 
of semiclassical quantization itself, and even yields numerically more accurate results than 
previous applications of filter-diagonalization (FD). The harmonic inversion method can be 
applied to the periodic orbit quantization of systems with both chaotic and regular classical 
dynamics, when the periodic orbit signal is calculated with Gutzwiller's trace formula []1],[2| 
for isolated orbits and the Berry- Tabor formula |§ for orbits on rational tori, respectively 
More generally, any signal given as a sum of 5 functions can be filtered analytically and 
analyzed using the methods described in Sections |II B| and |II Q For example, the technique 
can also be applied to the harmonic inversion of the density of states g(E) = J2 n $(E — E n ) of 
quantum spectra to extract information about the underlying classical dynamics fll2|j25| , |26l . 

The second method is the direct application of the Pade approximant (PA) to periodic 
orbit sums, and allows the resummation of the typically exponentially divergent terms of the 
semiclassical trace formulae. The Pade approximant can be applied when the total periodic 
orbit sum can be divided into partial sums with respect to an integer ordering parameter n, 
which can be related to, e.g., the cycle length of a symbolic code or the Maslov index of the 
orbit. 

The two methods for periodic orbit quantization have been demonstrated for systems 
with completely different classical dynamics, viz. the classically chaotic three-disk scattering 
problem and the integrable circle billiard. A detailed comparison of various semiclassical 
quantization methods for the three-disk repeller reveals that quantization by harmonic in- 
version and the Pade approximant to periodic orbit sums can even be applied in energy 
regions where the cycle expansion of the Gutzwiller-Voros zeta function does not converge. 

Evidently, the scope of applications of the methods presented in this paper is wide, 
and runs from simple models to challenging atomic, molecular, and mesoscopic systems 
with typically mixed regular-chaotic classical dynamics. Thus, even thirty years after the 
renaissance of semiclassics by Martin Gutzwiller, new and unexpected insights may be gained 
into the no man's land between classical and quantum physics. 
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FIG. 1. (a) Periodic orbit recurrence signal for the three-disk scattering system with R = 1, 
d = 6 without filtering. The signal in the region L < 35 consists of 93 non-equidistant periodic 
orbit contributions (including multiple repetitions), (b) Same as (a) filtered with frequency window 
w £ [182,218]. The band-limited signal consists of 201 equidistant data points with AL = 0.175. 
The solid and dashed lines are the real and imaginary part of C(L), respectively, (c) Semiclassical 
resonances obtained by harmonic inversion of the band-limited signal C(L) in (b). The dotted 
lines mark the borders of the frequency window. 
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FIG. 2. Semiclassical resonances (A\ subspace) for the three-disk scattering system with 
R = 1, d = 6. Squares: Harmonic inversion of the semiclassical recurrence signal; Plus symbols: 
Pade approximant to the periodic orbit sum; Crosses: Cycle expansion of the Gutzwiller-Voros 
zeta function pOj] . The dotted and dashed lines mark the borderline for absolute conver- 
gence of Gutzwiller's trace formula (Im k = —0.121557) and the Gutzwiller-Voros zeta function 
(Im k = —0.699 110), respectively. The Pade approximant and harmonic inversion results are found 
to converge deeper in the complex plane than do the Gutzwiller-Voros zeta function results. 



22 



OJ 



OJ 



'- (a) k=1 50-0.1 i 



10 
1 

0.1 
0.01 
0.001 
0.0001 
1e-05 
1e-06 
1e-07 

1e+08 
1e+06 
10000 
100 

1 - 
0.01 
0.0001 
1e-06 




H h-, 



-+— -+ 



1 1 r 

(b) k=150-0.5i 



i 7+ 








8 
n 



10 12 14 16 



FIG. 3. Convergence vs. exponential divergence of the partial periodic orbit sums for the 
three-disk scattering system with R = 1, d = 6 as functions of the order n at (a) complex wave 
number k = 150 — O.li and (b) k = 150 — 0.5i. Dashed lines and plus symbols: Error values 
e = \ A n — A n -i\ for the sequence A n without Pade approximation. Solid lines and squares: Error 
values £ = |^4^ A — -A^il for the Pade approximant to the periodic orbit sums. 
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FIG. 4. Some examples of periodic orbits of the circle billiard. The orbits are labeled by the 
numbers (m r ,m ip ) which correspond to the number of sides of the polygons and the number of 
turns around the center. The angle 7 is given by 7 = i:m v /m r . 
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FIG. 5. Real part (solid line) and imaginary part (dashed line) of the function l/g(k) for the 
circle billiard with radius R = 1 obtained by Pade approximant to the periodic orbit sum. The 
zeros agree perfectly with the exact positions of the semiclassical eigenvalues (from Eq. 36) marked 
by the squares. 
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